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Abstract 

In large extra dimensional Kaluza-Klein (KK) scenario, where the usual Standard Model (SM) 
matter is confined to a 3 + 1-dimensional hypersurface called the 3-brane and gravity can propagate 
to the bulk (D = 4 + d, d being the number of extra spatial dimensions), the light graviton KK 
modes can be produced inside the supernova core due to the usual nucleon-nucleon bremstrahlung, 
electron-positron and photon-photon annihilations. This photon inside the supernova becomes 
plasmon due to the plasma effect. In this paper, we study the energy- loss rate of SN 1987A due to 
the KK gravitons produced from the plasmon-plasmon annihilation. We find that the SN 1987A 
cooling rate leads to the conservative bound > 22.9 TeV and 1.38 TeV for the case of two and 
three space-like extra dimensions. 
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I. INTRODUCTION 



Recently it has been noted that the scale of quantum gravity, the four dimensional Planck 
scale Mpi(~ 10^^ TeV), is just a conjecture without much experimental support and the only 
experimentally verified scale of gauge interactions in four dimensions lies within the TeV 
scale. Therefore, the assumptions that gravitation becomes strong at the TeV scale, while 
the standard gauge interactions remain confined to the four dimensional spacetime, does 
not conflict with the today's experimental data These ideas solve the hierarchy problem 
without relying on supersymmetry or technicolour and the observed weakness of gravity at 
long distances is due to the presence of d new spatial dimensions large compared to the 
electroweak scale. This can be inferred from the relation between the Planck scales of the 
D = A + d dimensional theory Md and the four dimensional theory Mpi, which, for the 
toroidal compactification, is given by 

Mil = {2TxRfM'}j+\ (1) 

where R is the size of the extra dimensions. Putting Mn ~ 1 TeV then yields 

R ^ lof -^^cm. (2) 

For = 1, i? ~ 10^^ cm, this case is obviously excluded since it would modify Newtonian 
gravitation at solar-system distances. For d = 2, we get i? ~ 1 mm, which is precisely 
the distance where our present experimental measurement of gravitational strength stops. 
Clearly, while the gravitational force has not been directly measured beneath a millimeter, 
the success of the SM up to ~ 100 GeV implies that the SM fields can not feel these extra 
large dimensions, that is they are confined to only "3-brane", in the higher dimensional 
spacetime called "bulk" . Summarizing, in this framework the universe is D = 4 + c? dimen- 
sional with Planck scale near the weak scale, with d >2 new sub-millimeter sized dimensions 
where gravity and perhaps other fields can freely propagate, whereas the SM particles are 
localised on a 3-brane in this higher-dimensional spacetime. 

This theory predicts a variety of novel signals which can be tested using table-top experi- 
ments, collider experiments, astrophysical or cosmological observations. It has been pointed 



out that one of the strongest bounds on this physics comes from SN 1987A 
authors have done calculations to place such constraints on the extra dimensions 
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which we briefly discuss here. The graviton emission from plasmon-plasmon (photon in- 
side plasma of the supernovae (SN) becomes massive and is called as plasmon) annihilation 
might have deep impact on the supernovae cooling and can significantly alter the bounds 
on M/). Here we have investigated this possibility. This would be similar to Farzan's treat- 
ment of the Majoron emission in the supernova cooling process as a source of the upper 
bound on neutrino-Majoron coupling [3] and Raffelt's treatment on axion emission in pho- 
ton photon collision jsl. Several other mechanism for the SN 1987A cooling comprising the 
New Physics (beyond the Standard Model Physics) are already available in the literatute. 
Recently Das [9j and others (see [9] for other works) have explored the unparticle physics as 
a possible cooling mechanism of the supernovae SN 1987A, in which an unparticle stuff can 
be produced in the nucleon-nucleon bremstrahlung, electron-positron and photon-photon 
annihilations and thus cools down the temparature of SN 1987A. 



II. SUPERNOVA EXPLOSION AND COOLING 



Supernovae come in two main observational varieties: Type II are those whose optical 
spectra exhibit Hydrogen lines and have less sharp peaks at maxima (of 1 billion solar 
luminosities), whereas the optical spectra for the Type I supernovae does not have any 
Hydrogen lines and it exhibits sharp maxima [l^. Physically, there are two fundamental 
types of supernovae, based on what mechanism powers them: the thermonuclear supernovae 
and the core-collapse ones. Only supernovae la are thermonuclear type and the rest are 
formed by core-collapse of a massive star. The core-collapse supernovae are the class of 
explosions which mark the evolutionary end of massive stars (M > 8Mq). The kinetic 
energy of the explosion carries about 1% of the liberated gravitational binding energy of 
about 3x 10^^ ergs and the remaining 99% going into neutrinos. This powerful and detectable 
neutrino burst is the main astro-particle interest of core-collapse supernovae. 

In the case of SN 1987A, about 10^^ ergs of gravitational binding energy was released in 
few seconds and the neutrino fluxes were measured by Kamiokande ll| and 1MB [12] collab- 
orations. Numerical neutrino light curves can be compared with the SN 1987A data where 
the measured energies are found to be "too low" . For example, the numerical simulation in 
Q yields time-integrated values {E^J ^ 13 MeV, {EpJ ^ 16 MeV, and {E^J ^ 23 MeV. 
On the other hand, the data imply {EpJ = 7.5 MeV at Kamiokande and 11.1 MeV at 
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1MB [14j. Even the 95% confidence range for Kamiokande implies (-E'p^) < 12 MeV. Flavor 
oscillations would increase the expected energies and thus enhance the discrepancy ij]. It 
has remained unclear if these and other anomalies of the SN 1987A neutrino signal should 
be blamed on small-number statistics, or point to a serious problem with the SN models or 
the detectors, or is there a new physics happening in supernovae? 

Since we have these measurements already at our disposal, now if we propose some 
novel channel through which the core of the supernova can lose energy, the luminosity in 
this channel should be low enough to preserve the agreement of neutrino observations with 
theory. That is, Cnewchannei < 10^^ ergs s~^. This idea was earlier used to put the strongest 
experimental upper bounds on the axion mass |l5|. Here, we will consider the gravitons 
which can carry the energy from the core of the supernovae and escape into the bulk of 
the larger dimensional space. The constraint on luminosity of this process can be converted 
into a bound on the 4+d dimensional Planck scale Mr). Any mechanism which leads to 
significant energy-loss from the supernovae core immediately after bounce will produce a 
very different neutrino-pulse shape, and so will destroy this agreement, which in the case 
of a.io„ is exphcitly showu by Burrows's al fl. Raffe.t has proposed a s,mp.e analytic 
criterion based on detailed supernova simulations [19|]: if any energy- loss mechanism has 
an emissivity greater than 10^^ ergs g~^ s~^ then it will remove sufficient energy from the 
explosion to invalidate the current understanding of Type-II supernovae's neutrino signal. 
Similar arguments can be applied to other particles. The hypothetical majorons are one 
case in point |20|]. 



III. CONSTRAINTS ON EXTRA DIMENSIONS 



The most restrictive limits on come from SN 1987A energy-loss argument. If large 
extra dimensions exist, the usual four dimensional graviton is complemented by a tower of 
Kaluza-Klein (KK) states, corresponding to new phase space in the bulk. The KK gravitons 
interact with the strength of ordinary gravitons and thus are not trapped in the supernovae 
core. During the first few seconds after collapse, the core contains neutrons, protons, elec- 
trons, neutrinos and thermal photons(plasmons). There are a number of processes in which 
higher-dimensional gravitons can be produced. For the conditions that pertain in the core 
at this time (temperatures T ~ 30 — 70 MeV, densities p ~ (3 — 10) x 10^^ g cm~^), the 
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relevant processes are shown below 

• Graviton(^) production in Nucleon-Nucleon Brehmstralilung: N + N^N + N + Q 

• Graviton production in photon fusion: 7 + 7 ^ ^ 

• Graviton production in electron-positron annihilation process: e~e+ — > ^ 

In the supernovae, nucleon and photon abundances are comparable (actually nucleons 
are somewhat more abundant). In the following we present the bounds derived by various 
authors using nucleon-nucleon bremhmstrahlung and in the next section we give detailed 
calculation for photon-photon annihilation (including the plasma effect inside supernovae) 
to KK graviton process. We believe that although the dominant contribution will still follow 
from nucleon-nucleon bremsstrahlung, however, because of the large uncertainties involved 
in such a process calculation inside the hot plasma, the reliable bound will follow from 
plasmon -|- plasmon —>■ KK graviton process. It is worthwhile to mention here that in this 
work we have not considered the effect of plasmon width in the final continuum KK state 
production, which we believe if be taken into account will not substantiably change our 
bound on M^. We will not discuss the electron-positron annihilation to KK graviton as it 
does not give any significant bounds. 



A. Nucleon-Nucleon Brehmstrahlung 

This is the dominant process relevant for the SN 1987A where the temperature is com- 
parable to m^r and so the strong interaction between N's is unsuppressed. This process can 
be represented as 

N + N ^ N + N + g, (3) 

where N can be a neutron or a proton and ^ is a higher- dimensional graviton. 

The main uncertainty comes from the lack of precise knowledge of temperatures in the 
core: values quoted in the literature range from 30 MeV to 70 MeV. For T = 30 MeV and 
p = 3 X 10^^ g cm~'^, we list the results of various authors. 
CuUen and Perelstein ^ 

Md ^ 50 TeV, d = 2; (4) 
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Md > 4 TeV, d = 3; (5) 

Md > 1 TeV, d = 4. (6) 



Barger, Han, Kao and Zhang {4] 



Hannestad and Raffelt 



Md Z 51 TeV, d = 2; (7) 

Md ^ 3.6 TeV, d = 3. (8) 

Md > 84 TeV, d = 2; (9) 

Mb > 7 TeV, d = 3. (10) 



IV. METHODOLOGY OF CALCULATION 



Each KK graviton state couples to the SM field with the 4-dimensional gravitational 
strength according to [21 1 

^ = -^T. j d'x , (11) 



where the summation is over all KK states labeled by the level n. Here k = ^/T6^TGN and 
Gn = ^/Mpi, the 4-dimensional Newton's constant. T^^^ is the energy-momentum tensor of 
the SM and h^"'^'^ the KK state. 

Since for large R the KK gravitons are very light (because ~ they may be 

copiously produced in high energy processes. For real emission of the KK gravitons from 
the collision of SM fields, the total cross-section can be written as 



CTtot = «:'E^(^) ' (12) 



n 



where the dependence on the gravitational coupling is factored out. Because the mass 
separation of adjacent KK states, 0{1/R), is usually much smaller than typical energies in 
a physical process, we can approximate the summation by an integration according to 

E - / P{Od{mll (13) 

where the density of KK states p{mV) = ^+tf 4d^3d/2p(^/2) (^l)''^"^''^^- Here we have used the 
relation M|,; = (27ri?)°'M|,+'^. 



Now for a generic 2 N body scattering, the scattering cross section is given by 

where Flux = 4EiE2Vrei- Here Ei, E2 are the energies of the initial particles 1 and 2 whose 
masses are mi and m2, respectively and Vrei is the relative velocity between them. 
For a general reaction of the kind a + 6 — > c, the above expression takes the form 

a=^\M^''^^KS~rnl). (15) 

In the center-of-mass frame, we use the notation for the total initial energy, 

^fS = Ei + E2 (16) 

FluX = AEiE2Vrel = ^^\^. (17) 

where |p| = |pi| = |p2| = ^ ' ^2'%/^'"^'^'' ^"^^ are the energies of the particles a and 

b. The function A(x, y, z){= + y'^ + z'^ — 2xy — 2yz — 2zx), is the standard ii'cillen function. 

Since we are concerned with the energy loss to gravitons escaping into the extra dimen- 
sions, it is convenient and standard Q to define the quantities ia+b^c which are the 
rate at which energy is lost to gravitons via the process a + b c where c has a decay width, 
per unit time per unit mass of the stellar object. In terms of the cross-section aa+b^c the 
number densities na,b for a,b and the mass density p, e is given by 

_ {nanb(y(^a+b-,c)VrelEcm) 
^a+b~^c. — UoJ 

where the brackets indicate thermal averaging and Ecm{= Ea + Eb) is the center-of- 
mass(c.o.m) energy of the two colliding particles a and b. Note that in the present case 
the final state KK graviton, although has smaller decay width but is stable over the size of 
the neutron star because of it's large life time ~ 10^(100 MeV/m)^ yr (See [21I) and thus it 
can escape the supernovae while allowing it to cool. 

V. GRAVITON PRODUCTION IN PLASMON FUSION 

Photons are quite abundant in supernovae. Due to plasma effect inside the supernovae, 
photons becomes effectively massive. These massive photons(of mass mA,say) are known as 
plasmons. Our interest is in the plasmon-plasmon annihilation to KK graviton i.e. 

'jp{ki)+-fp{k2)^KK{p). (19) 



The plasmon-plasmon-graviton {G'^^{q)A'^{ki)A'^ "^(^2)) vertex [2l| is given by 



X, 



IK 
' 2" 



[m\ + ki.k2)Ci 



(20) 



where the symbols Cpu,pa, Dp,y^p„{ki, ^2) , E^y^p„{ki, ^2) are defined as 



fj,u,pcr 



Dpu,pa{ki, h) = Vfiuhahp - 'np„kiyk2p + Vpphahu - Vpahf^hu + [fi ^ Jy) 
Epu,p^iki, ^2) = Vf^uikipha + hpha + hpha) - [riv<ykif,kip + riypk2p,k2a + ^ ^)] ■ 
Here we work in the unitary gauge(.^ 00). In the c.o.m frame, the momentum vectors for 
this reactions are 

(21) 
(22) 

= {Eg, 0,0,0). (23) 
It often turns out to be more convenient to keep the polarizations explicitly. The polarization 



fcf = (^i,0,0,fc), 
{E2,0,0,-k), 



^2 



vectors 



2l[ | of a massive graviton are 

- 2e^6± 



p 5 



V ~ V2(e^e^ + e^e^j , 



-^pu 



3 [^p 



Here and e° are the transverse and longitudinal polarization vectors of a massive gauge 
boson. For a massive vector boson(e.g'. plasmon) with momentum k^ = {E, 0, 0, k) and mass 
rriA, 

1 



e+(/c) = ^{0,1,1,0) 
^~p{k) = ^(0,-1, z,0) , 
iik) = —{k,0,0,-E) . 



(24) 
(25) 
(26) 



The plasmon and graviton polarization vectors satisfy the following normalization and po- 
larization sum conditions 



^S-', j:e;{k)et*{k) = -vp.+ ^'^'' 



m 



2 ' 



(27) 
(28) 



s=l 



where Bf,yp^{j>) is given by 

Pt^Pp \ ( PuPa 



Bf^upAp) = 2 r/^p - 



2 2/1 i'^'^ 1 



n n , 



The total squared amphtude, averaged over the initial three polarizations (since massive 
plasmons have three state of polarizations) and summed over final states for the process 
7p(^i) + 7p(^2) Gkk{p) is 

2 ,^2 



l-^r = (^j Y.\^f = Y2 + + ^3 + T| + Tt) , (30) 

where {i = 1, ..5) are given in appendix A. Substituting this in ( |T5l) and using ( |T6l) and 
f fT7|) . the total cross-section cry for this process is obtained as 

1 /■ 

(Tt = 5I^7P7P-G,,(5,m^) = P("^^) ^iS-ml) \Mf d{ml) 

n 

~ 9 4'^7r-r(rf/2) 1^M2 J ^ ^ 

where ^ = -1 + f and AT = ^ (^|^ + + f + 16^ + f ). While deriving 

Eq. [311 "W'e have used pinifi) = (^4^^Y/'^T{d/2) Planck scale relation Eq. [TJ 

The volume emissivity of a supernova with a temperature T through this process is 
obtained by thermal-averaging over the Bose-Einstein distribution. Hence, the energy loss 
rate (e^p = '^~Q'yp) due to plasmon plasmon annihilation is given by (similar to that of the 
energy loss rate via 77 —>■ vv 



1 1 



PSN TT Ju)o 



doji 



pUJi/T _ 1 



duj- 



2U/2 



S{uJi + UJ2] 

2lUiUJ2 



(TT, 



(32) 



/ 2 2\l/2 

where ax is given in Eq. [HU Note that N^p = du '^^'^^^^^'^'^ — is the number density 
of thermal photons, or rather of transverse plasmons. In the present case, we treat the 
plasmon to be transverse (with the dispersion relation given by c<j^ = uJq + |k since the 
contribution coming from the longitudinal plasmon is typically smaller 23|, 
above, Uq corresponds to plasma frequency in the supernovae core. 



25|. Also in 
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Finally introducing the dimensionless variables Xi = uJi/T{i = 0, 1, 2) and taking (the 
transverse plasmon mass) to be equal to uq, we rewrite the above Equation as 

^7P - , ^2+d_A / :^T7t / —7^ — i-idj 



P57vM|+V7xo 7x0 xiX2 

where 

1 1 r 4 T2 2 16m^ 1 16m^ 1 17 
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VI. NUMERICAL ANALYSIS 

The SN 1987A energy loss due to KK graviton emission produced in massless photon- 
photon annihilation already put some bound on the effective scale of gravity Md for d = 2 
and 3 (see j^). Here we study the modification of the above bound in a scenario where the 
plasma effect on photon is taken into account. In our analysis, the key working formula is the 
Eq. [33] which describes the supernovae energy loss rate due to plasmon(7p) + plasmon(7p) 

KK graviton(G'7^i^). Now for any kind of cooling mechanism which corresponds to an 
emissitivity > 10^^ erg g^^ s^^ would invalidate our current understanding of Type-IIA 
supernovae's neutrino signal. So the consistency with the neutrino signal requires the energy 
loss rate < 10^^ erg g~^s^^. This gives rise the lower bound on M/). In Fig. 1 we have shown 
the energy loss rate to KK gravitons as a function of the scale Mr, for different number of 
extra dimensions d. The right and left curves respectively stands for d = 2 and 3. In this 
plot, the inputs taken are as follows: ujq = m^(plasmon mass) = 19 MeV, the supernovae 



temperature T = 30 MeV and the supernovae core density p ^ 10^^ g cm^^ 19|]. The 
horizontal line corresponds to the upper bound on the supernovae energy loss rate. The 
intersection of this curve with the other two gives rise the following lower bound on M^: 
for = 2 we find > 22.9 TeV, whereas for d = 3, Md > 1.38 TeV. The bound on Mn 
as obtained here is somewhat stronger((i = 2) and weaker((i = 3) than that obtained in ^ 
which are 15 TeV and 1.6 TeV, respectively for d = 2 and 3, where the relevant process of 
interest was photon-photon annihilation to KK gravitons. Also note that the bound on M^, 
that we find from plasmon-plasmon annihilation to gravitons is somewhat weaker than the 
one obtained from the nucleon-nucleon brehmstrahlung which are 51(3.6) TeV for = 2(3), 
respectively. Finally, note that the present supernovae SN 1987A cooling analysis does not 
allow us to put any bound on Md for c/ > 4. 
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Fig. 1. The supernovae energy loss rate dt/dt (erg g^^s^^) due to KK graviton emission 
produced in plasmon plasmon annihilation is shown as a function of M£,(GeV) in Fig. 1. 
For the right curve d = 2, whereas for the left d = 3. The upper horizontal curve corresponds 
to de/dt < 10^^ erg g~^s-^. 

VII. CONCLUSIONS 

In summary, we found that the emission of KK graviton by plasmon-plasmon annihilation 
from SN 1987A puts the conservative bound on the effective scale M/? of the large extra 
dimensional model in the case oi d = 2 and 3. Taking a conservative estimate of the 
supernovae temperature T = 30 MeV and plasmon mass = 19 MeV (equal to the core 
plasma frequency cUo), we find Md > 22.9 TeV for d = 2 and Md > 1.38 TeV for d = 3. No 
bound on follows from the present analysis for d > A. 
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3 V mi mt^ 
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